On the non-V-representability of currents in time-dependent many particle systems 
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We argue that an arbitrarily chosen time-dependent current density is generically non-V- 
representable in a many-particle system, i.e., it cannot be obtained by applying only a time- 
dependent scalar potential to the system. Furthermore we show by a concrete example that even 
a current that is V-representable in an interacting many-particle system may (and in general will) 
turn out to be non- V-representable when the interaction between the particles is turned off. 
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I. INTRODUCTION 

Since the beginning of density functional theory 
(DFT>ii2i^, the problem of answering the following ques- 
tions has been recognized to be of fundamental impor- 
tance: 

(1) Given a (positive) particle density n{r) in an N- 
particle system is there a local potential V{r) that 
produces this density in the ground-state of the sys- 
tem? 

(2) If a certain particle density n{r) arises in the 
ground-state of an N-particle system subjected to a 
local potential V{r), is there a local potential Vs{r) 
that produces the same density in the ground-state 
of the same system with the particle-particle inter- 
actions turned off? 

The first question is known as the V-representability 
question and a given density is said to be V-representable 
if the answer is affirmative. The original formulation 
of DFT by Hohenberg and Kohni made heavy use of 
the assumption that "all reasonable densities are V- 
representable" . Subsequent worki has shown that this 
assumption is necessary only when one tries to prove the 
existence of the functional derivative of the energy func- 
tional. The second question lies at the very heart of the 
Kohn-Sham formalism:- Recall that within this formal- 
ism one tries to obtain the ground-state density of an 
interacting many-particle system by applying a local po- 
tential to a non-interacting version of the same system. 
Obviously, it is vital to the success of this theory that the 
target density, which is by assumption V-representable, 
be also non- interacting V-representable - i.e., the answer 
to question (2) must be affirmative. 

Even though mathematically rigorous answers to the 
two V-representability questions are not known, DFT 
has been widely applied to the calculation of the elec- 
tronic structure of matter. In these calculations it is 
tacitly assumed that the set of V-representable densi- 
ties in both interacting and non-interacting systems is 
dense enough to approximate to an arbitrary level of ac- 
curacy any physical ground-state density (these beliefs 
are supported by mathematical work on lattice system^) . 
These assumptions have been automatically transplanted 



to the relatively younger field of time-dependent density 
functional theory (TDDFT)SiL&, where the questions are 
whether a given time- dependent particle density n(r, t) 
evolving from a given initial state 4* can be produced by 
a local time-dependent scalar potential ^(r, t) and, in 
the affirmative case, whether the same density can also 
be produced by a time-dependent scalar potential Vg (r, t) 
starting from an initial state now with interactions 
turned off. 

In this paper we are not going to challenge the wis- 
dom of the standard V-representability assumptions in 
DFT or TDDFT, but rather examine whether such as- 
sumptions can be plausibly extended to the particle cur- 
rent density j (r , t) . There are good reasons to under- 
take this study. During the past ten years we have 
seen many indications that the time-dependent current 
densityjSiifliiiiiSii^ together with the initial state of the 
system (and hence the initial density) may provide a 
more fundamental description of the dynamics. Indeed, 
the proof that the time-dependent current density deter- 
mines the external scalar potential is the very first step in 
the proof of the Runge-Gross (RG) theorem - the foun- 
dation theorem of TDDFT. ^ However, the RG theorem 
does not say anything about the V-representability ques- 
tion for the current, i.e. whether a given time-dependent 
current can be produced by a local time-dependent scalar 
potential. Reasoning by analogy with the particle den- 
sity this has led (or, as we are going to show, misled) 
some workers to believe that the V-representability as- 
sumption for the current density is about as plausible as 
the corresponding assumption for the density, and that 
therefore any physical current density can be approxi- 
mated to an arbitrary degree of precision by the current 
density generated by suitably chosen scalar potential, in 
an interacting as well as in a non-interacting system^ 

The purpose of this paper is to show that this is not 
the case. Due to the vector character of the current den- 
sity it is usually impossible for an arbitrary current to 
be generated by a single scalar function - the potential. 
Even in those special (but physically very relevant) cases 
in which this can be done, the current will not be simul- 
taneously representable in the non-interacting system. A 
more general theory that makes use of an effective vec- 
tor potential to generate the current is therefore needed: 
such a theory exists and it is known as time-dependent 
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current density functional theory (TDCDFT)2iH. 

II. NON-V-REPRESENTABILITY OF GENERIC 
CURRENT DENSITIES 

In this section we present our main argument against 
V-representabihty of the current density. Recall that the 
vector field j(r, t), like any other vector field, can be writ- 
ten as the sum of a longitudinal component jL{f,t)^ a 
transversal component Jt('^j^)i and a constant that can 
be assumed to vanish if the full current has to vanish at 
infinity: 

'][r.t)^]L{r,t)+Mr.t). (1) 

The longitudinal current density {r, t) is curl- free and 
can therefore be represented as the gradient of a scalar 
field, while the transversal current jT{f,t) is divergence- 
free and can therefore be represented as the curl of a vec- 
tor field. The spatial Fourier components of jl and jt, 
denoted by ^^((f, t) and jT('f, t) respectively, are obtained 
by projecting the Fourier component of j{q,t) along di- 
rections parallel and perpendicular to the unit vector q, 

jLiq,t) ^ {j{q,t) ■ g, (2) 

and 

jriq, t) = i{q, t) - jL{q, t). (3) 

Notice that the particle density is completely determined 
by the longitudinal current density since, according to the 
continuity equation, one has 

^ = -V.J.(f,t) (4) 

with initial condition n{f,t) — noif). 

We begin our argument by assuming that a certain 
current density j{r, t) is V-representable, and let jL(r, t) 
and Jt [r, t) denote its longitudinal and transverse com- 
ponents, respectively. According to the RG theorem, the 
potential V{f^t) that produces j{f,t) is unique up to 
an arbitrary function of time. Consider now a second 
current density = j'^{r,t)+j^{r,t) = jL{r,t) + 

i'rp{r^t), i.e. j'ir,t) differs from j{f,t) only because its 
transverse component, j'rp{f,t) differs from jT{r,t). We 
claim that j'{r,t) is not V-representable. Indeed, if it 
were V-representable, then there would be a potential 
V'{r,t) 7^ V{r,t) that generates it (here and in the fol- 
lowing the 7^ sign means that two potentials differ by 
more than a mere function of time). But this is im- 
possible, since, according to Eq. Q), these two different 
potentials would give the same particle density, in con- 
tradiction with the Runge-Gross theorem.^ Thus, for a 
given longitudinal current density {r, t) there is at most 




FIG. 1: A simple cartoon of the sparsity of the V- 
representable current densities. The V-representable current 
density j lies on a continuous hypersurface (here schematized 
as a continuous curve) in current density space. Due to the 
continuity of the mapping between jl and jr a sufficiently 
small neighborhood of the non- V-representable current den- 
sity j' contains only non- V-representable current densities. 

one transverse current density jx {r, t) such that the full 
current density jhir, t) + jrir^ t) is V-representable. 

The ease with which, given a V-representable cur- 
rent density, we were able to construct infinitely many 
non- V-representable ones is a strong indication that V- 
representable currents are a rather exceptional occur- 
rence in the space of all possible currents. To strengthen 
the argument let us make the plausible assumption that 
the mapping from potentials to V-representable currents, 
via the solution of the time-dependent Schrodinger equa- 
tion, is not only invertible (RG theorem), but also con- 
tinuous. This implies that, within the subset of V- 
representable current densities, jTir,t) is a continuous 
functional of ji(f, i). Let j{r,t) ~ iL{r,t) + jTir,t) be 
a V-representable current density. Consider then a small 
"neighborhood" of the non- V-representable current den- 
sity j'(r, t) = jL{r, t) -t- %{r, t) and let Ji(r, t) be a (hy- 
pothetical) V-representable current density in this neigh- 
borhood. Since, by choice, the longitudinal component 
of ii, JiL, is close to jL, the continuity of the mapping 
from to jT for V-representable currents implies that 
the transverse component of ji, jit, is close to jt- But, 
this cannot be true for a sufficiently small neighborhood 
of j'(r, t) if the difference between jt and is finite (see 
Fig. 1). We conclude that every non- V-representable 
current density is surrounded by a neighborhood that 
contains only non- V-representable current densities: the 
set of V-representable current densities is not "dense" in 
the space of all possible current densities. In Appendix IXI 
we give another more mathematical argument to support 
our statement. 



III. INTERACTING VERSUS 
NON-INTERACTING V-REPRESENTABILITY 

Undeterred by the above arguments one might insist 
that, after all, the task of the Kohn-Sham theory is to 
approximate V-representable current densities by non- 
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interacting V-representable ones. We know that the set 
of V-representable current densities is characterized by a 
certain jx vs Jl relation: it is the presence of this con- 
straint that makes the set so "sparse" . Similarly, in the 
non-interacting system, the set of the V-representable 
current densities is characterized by another jx vs Jl 
relation. Wouldn't it be nice if these two relations hap- 
pened to be the same relation, so that a Kohn-Sham 
potential yielding the correct density would also auto- 
matically yield the correct current density? 

This conjecture has found its way in the literature^, 
so it is important to examine it carefully. In this sec- 
tion we construct an example of a V-representable cur- 
rent density, which is definitely non- V-representable in 
the non-interacting (Kohn-Sham) system. Thus, in this 
example the jx vs Jl relation of the interacting system 
turns out to be incompatible with jx vs jx relation of 
the corresponding non-interacting one. We will argue 
(without proof) that this state of affairs is quite generic 
for currents generated by scalar potentials in anisotropic 
interacting systems. 

The system we consider is a two-dimensional interact- 
ing electron liquid which is initially in the ground-state 
with inhomogeneous density 

no{r) = h{l + 2jcosG ■ r), (5) 

where G is a two-dimensional vector parallel to the x- 
axis, n is the average density, and hj the amplitude of 
the density modulation, with 7^1. This density is 
produced by the application of the static potential 

Vo{f)^^^cosG-f (6) 

to an initially homogeneous electron liquid of density ni^ 
Here x''(G) is the static density susceptibility of the ho- 
mogeneous electron gas of density n at wave vector Gi^ 
We now apply to this system a time-dependent, periodic 
scalar potential of the form 

Vif,t) = t/e^^'?-^--'^*) + c.c. , (7) 

where q is a two-dimensional wave vector, which we as- 
sume to be much smaller in magnitude than both G and 



kp {kp being the Fermi wave vector corresponding to 
the average density n). The frequency lo is assumed to 
be larger than both vpq and vpG, where vp is the Fermi 
velocity associated with kp. 

It is well known that the time-dependent potential of 
Eq. |7J) can be recast as a longitudinal vector potential 

A[r,t) = i^e*(9-"--*) + c.c. , (8) 

so we can say that the induced current density is 

j„(f,t) =^XaM9 + ^G,g,c^)*^e*[('"'+'<3)---*l + c.c. , 

1,13 

. ^ . . . . 

where Xaf3{<l + ^G, q, oj) is the dynamical linear response 
function of the inhomogeneous liquid, which connects the 
cartesian (3 component of the vector potential amplitude 
at wave vector q to the cartesian a component of the 
current density amplitude at wave vector q -\- IG, where I 
is an integer. It should be noted that this response func- 
tion coincides with the homogeneous response function 
xj^^((f, w) up to corrections of order 7^ when I = 0, and 
is of first order in 7 when ^ = ±1. The components with 
^1 > 2 are of order 7^ at least. 

We now want to show that the exact current of Eq. @ 
cannot be obtained in a non-interacting system subjected 
only to scalar potentials that yield the exact density. We 
first notice that the ground-state density hq (r) is enforced 
in a non-interacting electron gas by the scalar potential 

where Xo(^) is the static density susceptibility of the 
noninteracting electron gas of density n at wave vector G. 
We know from the invertibility of the mapping between 
vector potentials and currents that the exact current den- 
sity of Eq. © can be generated in the non-interacting 
electron gas by one and only one time-dependent vector 
potential. 



AsAr,t) = A^{r,t)+ ^ /^^^ ^^(g + ;g, g + ^G, ^)X/35(g ?'G, c^)^e't(«+'^)-^^-"*l c.c, (11) 

where fHxc,ap{q + ^G, q + I'G, uj) is the sum of the so-called exchange-correlation kernel fxc,aii{q + IG, q + I'G, ui) and 
the Hartree kernel /jj,a/3(9 + IG, q + I'G, t^) = p^^^(g| '■'^"'^"^"^a'^'^^"^'^"' ^if for our system. Thus, in order to prove that 
the current of Eq. @ cannot be generated, in the non-interacting system, by a scalar potential, we only need to show 
that the vector potential As{r,t) has a finite transversal component. For, if this is the case, then the uniqueness of 
As guarantees that the current density cannot be produced by a purely longitudinal vector potential, and hence not 
by a simple scalar potential. 
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In order to establish the existence of a transversal component of Ag we focus on the Fourier component at wave 
vector q and we discard both A and the contribution of the Hartree kernel because they are purely longitudinal fields. 
The quantity of interest is thus the exchange-correlation vector potential 



— 

LI3,S 

The first non- vanishing term in the expansion of this quantity in powers of 7 is of order 7^: 

A^cAq^^) - ^{fL.ai}{q,^)x^psi^^?^^) + fxlapi?^?^^)x%{q^ 



(12) 



(13) 



where the quantities with superscripts and refer to the inhomogcncous sytcm of density no{r) and are first 
order and second order in 7, respectively. 

To proceed we now make use of certain exact identities, which can be obtained starting from identities that were 
derived in Ref. from a careful consideration of the behavior of the current response function and its associated xc 
kernel under transformation to an accelerated reference frame. These "acceleration identities" are summarized in the 
following four equations: 



xS(G,0,^) = -i 



xS(o,o,^) 
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x\G,uj) 


TO 


5af3 - 


xHG) 








y muj ^ 


1 xHG) 



Ptp{G) 



x\G,Lo) 



X''{G) 



PapiG), 



/ic!a/3(G, 0, u) = { [f^,^{G, u) ~ f^LiG, 0)] Ptp{G) + ftAG, io)Plp{G)] , 



/rc:a/3(o,o,^) = -27/£;„^(G,d,^), 



(14) 
(15) 

(16) 
(17) 



where f^cL^G, io) and /^ct(G, ^) are the longitudinal and 
transverse xc kernels of the homogeneous electron liquid 
at density n and wave vector G, while P^piG) = -^^frr^ 



and P^JG) 



PadiG) are the projectors paral- 



^p\^l = 5ap - -L ^py 

lei and perpendicular to the direction of G, respectively. 
The derivation of these identities is briefly presented in 
Appendix ^ 

What makes these identities relevant to the evalua- 
tion of Axc{qi^) is the fact that in our model q is much 
smaller than \G\ or kp: therefore the quantities appear- 
ing in Eq. (|13|l can be evaluated in the g — > limit where 
they reduce precisely to the quantities that appear in 
Eqs. I|14|) ~ (|17|l . Underlying the calculation is of course 
the assumption that the g — > limits of the current re- 
sponse functions and xc kernels are regular - an assump- 
tion we have presently no reason to doubt. 

By making use of Eqs. lC3I)-(C7)) in Eq. and by 

discarding all but the leading-order terms in q we arrive, 
after some algebra, at 



Axc{q,^^) 



2j 



^nG{G-q)x''{G,uj) 



x'^iG) 



X [/^,JG,c.)-/,^,i(G,0)] 



V 



(18) 



This vector has a component perpendicular to q, unless 
q happens to be either parallel or perpendicular to G: 
choosing q in any other direction provides the desired 
example of non V-representability. 

It is also interesting to ask what is the behavior of the 
Fourier components of A^c at wave vectors q±G. These 
components are first order in 7 and one might think that 
they lead more directly to the desired result. Remark- 
ably, this is not the case: A calculation very similar to 
the one described in the previous paragraphs reveals that 
these components are purely longitudinal, i.e., parallel to 
q±G. One needs to go to at least second order in 7 to 
see a transverse component of A^c- 



IV. DISCUSSION 

The example worked out in the previous section shows 
that a perfectly legitimate V-representable current den- 
sity can turn out to be non- V-representable in the non- 
interacting system. We believe that this state of affairs is 
generic. Only in exceptional cases will the current den- 
sity be V-representable in both the interacting and the 
non-interacting versions of the same system. Hence, in 
general, the Kohn-Sham equation does not give the cor- 
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rect value of the transverse current density. 

We may now ask, how big an error does one make if one 
insists on calculating the transverse current by means of 
the Kohn-Sham theory? Going back to our model system 
it is not difficult to see that the Fourier components of 
the Kohn-Sham potential are given by 



VKs{q + lG,uj) = 



\q + lG\ 



\A,^Liq + lG,io)\ 



(19) 



i.e., the Kohn-Sham potential is simply the scalar rep- 
resentation of the longitudinal part of the vector poten- 
tial As calculated in the previous section {As,l and As^t 
are the longitudinal and transversal component of Ag, 
respectively). The above equation is accurate up to cor- 
rections of order 7'^. Indeed, because Ag^x is of order 7^ 
its influence on the longitudinal current begins at order 
7^^, implying that alone fully accounts for the lon- 
gitudinal current density up to corrections of order 7'^. 
On the other hand, because Vks is equivalent to As^l it 
clearly fails to produce the part of the transverse current 
that is due to Ag.r- This is of order 7^ for the Fourier 
components at wave vector q + lG with I ^ and of order 
7^ for the Fourier component at wave vector q. We con- 
clude that the error on the transverse current is overall of 
order 7^: this may partly explain the difficulty of finding 
examples in which the Kohn-Sham current density differs 
significantly from the exact one. 

Where does this leave us as to the application of the 
time-dependent Kohn-Sham theory to the calculation of 
current densities? From a fundamental standpoint it is 
clear that only the time-dependent CDFT can provide 
the right answer. In time-dependent CDFT one does not 
need V-representability, but only the much weaker A- 
representability assumption. We know that this assump- 
tion holds true in linear response theory, and it is highly 
reasonable to assume that the set of A-representable cur- 
rent densities is dense in the space of all current densities. 
On the other hand, we have also found that the error en- 
tailed by the use of the ordinary Kohn-Sham equation of 
TDDFT is of second order in the parameter that mea- 
sures the strength of the density non- uniformity, and may 
perhaps be reasonably neglected in practical implemen- 
tations of the theory that are based on the local density 
approximation. 
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APPENDIX A: A FORMAL PROOF OF THE 
NON-DENSE V-REPRENSENTABLE CURRENT 
DENSITY SET 

In this section we discuss a different proof of the fact 
that the V-representable current density set is not dense 
in the set of physical currents. Let us recall that, math- 
ematically, a subset X of the set A is dense in A if and 
only if every element of A is the limit of a sequence of 
elements in Xi^lt is easy to see that ii X A and X 
is closed, i.e. any sequence of elements in X has limit in 
X, then X cannot be dense in A. 

Our first observation is that the physical currents, 
which are integrable and diffcrentiable together with 
their derivatives, form therefore a Sobolev space. Let 
us now construct a sequence {ji } of V-representable cur- 
rent densities and let us study its limit I —> 00. Be- 
cause the Sobolev space is closed, the limit joo{r,t) = 
limi^oc jii^,t) is still an integrable and differentiable 
fmiction, i.e., it is still a physical current. This is also true 
for the corresponding densities {ni} and wave-functions 
l^";}. We then consider the sequence of potentials {V;} 
that generates the sequence of current densities {ji}. The 
hypothesis of the continuity of the mapping between po- 
tentials and currents formulated in the text ensures the 
existence of the limit 



lim Vi{f,t)^Vao{r,t). 



(Al) 



We use this result in the sequence of Schrodinger equa- 
tions 



(A2) 



where Hi is defined by the potential Vi, to pass to the 
I —> 00 limit and arrive at the conclusion that is V- 
representable, and the potential which defines it is Voo- 
Thus all the sequence of V-representable currents have 
a V-representable limit, i.e. this set is closed. Because 
we proved that the set of non- V-representable currents is 
not empty, the V-representable current densities do not 
form a dense subset of all the physical current densities. 



APPENDIX B: ACCELERATION IDENTITIES 
FOR THE CURRENT RESPONSE FUNCTIONS 
OF A WEAKLY INHOMOGENEOUS ELECTRON 
LIQUID 

Our starting point is the acceleration identity in real 
space 



Ssp + 



V'.V^K.(r") 
m(iw)2 



dr = dafj, 



m 



_ . . (Bl) 
which was derived in Ref. 1141 beginning on page 205. 
Both sides of this identity can be expanded in a series 
of 7: The Fourier transform of the coefficients of this 
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expansion will give us the identities H14(l - H15() . The first 
order in 7 evaluated at wave vector G gives: 



/( vUJ lib 
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Inserting 



Vo(G) 



717 



(B2) 
(B3) 



from Eq.@, and recalling that 

Y,xUG.u;)GsG^ = x\G,u)Lo'P^^iG) , (B4) 



where x'*(G, w) is the density-density response function, 
we arrive at Eq. (|14|) . 

To derive Eq. p5|) we take the q — component of both 
sides of Eq. HB1|) of the second order in 7. Since Vo(r) 
and no{r) are given by Eqs. jSJl and © up to corrections 
of order 7^^ we readily obtain 



xS(o,o,w) 



-E[x«i(0'^'^) + ^«i(0,-G,< 
s 

GsGp 



muj 



-V„iG). 



(B5) 



Inserting Eq. (14) for x^"'^ (0, G, and Eq. (|B3p we im- 
mediately arrive at Eq. (|15|l . 



We proceed similarly for the last two identities, 
Eqs. H1B|I and p7ll . The starting point in this case is 
Eq. (27) of Ref. ,14: 



fxc,ap[r,r ,w)no(r )dr = ^ , (B6) 



where, to the required accuracy^^ 

V,M^ = 2njfi'^L{G, 0) cosG • r (B7) 

is the exchange-correlation part of the static Kohn-Sham 
potential K(^- Notice that /;^cl(^'0) is the scalar 
exchange-correlation kernel of the homogeneous electron 
liquid at zero frequency, quite different from the tenso- 
rial and frequency-dependent exchange-correlation kernel 
of the inhomogeneous system, which appears on the left 
hand side of Eq. I|B6(I . Taking the Fourier component at 
wave vector G of both sides of Eq. ljB6|) to first order in 
7, and recalling that 



fxc,af3{G,Uj) 



G^ r 

Cj2 



/^,^(G,c.)Pi^^(G) 



+ /,^^(G,a;)Pj^(G) 



(B8) 



we arrive at Eq. (|15|l . Taking the Fourier component of 
both sides of Eq. (|B6p at wave vector to second order 
in 7 we finally arrive at Eq. H17|l . 
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